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ABSTRACT 

Until now the existence of extended emission is an intriguing property of short bursts. 
It is not clear what is the nature of the extended emission. It might be a rising x-ray 
afterglow, or it could be a manifestation of the prolonged activity of a central engine. 
We consider short duration gamma-ray bursts, emphasizing the common properties 
of short bursts and short burst with extended emission. Assuming that the extended 
emission with broad dynamic range is a common property of short bursts, we propose 
a two jet model which can describe both short main episode of hard spectra emission, 
specific for short bursts, and softer spectra extended emission by different off axis 
position of observer. The toy model involves a short duration jet powered by heating 
due to vv annihilation and long-lived Blandford-Znajek jet with significantly narrow 
opening angle. Our proposed model is a plausible mechanism for short duration burst 
energization, and can explain both short burst with and without extended emission 
within a single class of progenitor. 
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1 INTRODUCTION 

It is commonly accepted that Gamma-Ray Bursts (GRB) 
consist of two populations - long and short duration bursts. 
The presence of separate short Gamma-Ray Burst (SGRB) 
population was suggested in (jMazets fc et al. 19811, and 



. . confirmed by the BATSE experiment (Kouveliotou & et al. 
> fl993b . The nature of S GRB can be merging of compact com- 
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panions in close binary systems such as neutron stars (NS) 
or N S Black-Hole (BH) ( jBlinnikov eFaLl[l984} |Paczynski| 
|1986[ ). Short GRBs have several distinct phenomenological 
properties which we briefly discuss. 

SGRBs populate the short mode of a bi-modal dura- 
tion distribution. The duration parameter, Tgo, is defined 
as continuous interval comprising the 90% of GRB emission 



in gamma-ray domain ( Koshut fc et al.|1996 l. It is also ac- 
cepted that the Tgo < 2 s is a good dividing line for short 
and long bursts for mid-energy-range (30 — 300 keV) ex- 



periments e.g. ( Kouveliotou fc et al.|1993 l. Later on a more 
complex criterion has been suggested in order to distinguish 



phenomenologically short and long GRB (Donaghy & et al. 
2006||Zhang fc et al.|2009 l 



SGRBs have a harder spectrum than long bursts ( Kou- 
veliotou fc et al.|1993[ ), and there is no spectral lag in their 



light curves in comparison with long bursts, where light 
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curves of the same GRB in soft channel lag the light curve 
in harder channels ( |Norris et al.|2000 1 . 

In a new era of rapidly slewing robotic telescopes, the 
optical afterglows of SGRBs has been detected, and in some 
cases it is possible to determine a redshift, z, of some SGRB 
sources. It turn out that most SGRBs have z < 1 and equiv- 
alent isotropic energy Ei so in the range 10 48 to 10 51 ergs 
( Gehrels et al.|2009| . 

Finally, until now there has been no cases of a supernova 
signature detected, either spectroscopically, or photometri- 
cally. 

One of most intriguing properties of SGRBs is that they 
exhibit extended emission. Prompt gamma-ray emission of 
SGRB consists of a short main episode, sometimes resolved 
into substructure of short pulses, which we call the Initial 
Pulse Complex (IPC). The duration, T 90 , of the IPC of the 
short burst is usually less than ~ 2 s. However, in the sum 
of light curves of many short bursts, aligned relative to the 
main peak of the IPC, significant extended emission (EE) 
has been observed up to ~ 100 s in different experiments 
and energy bands, in particular the BATSE 25 — 110 keV 
( Lazzati et al.|2001[) and 50 - 300 keV ( |Connaughton|2002 \ 
(see Figure|l[) (|Connaughton|2002[ ), Beppo-Sax 40-700 keV 
(jMontanari et al.||2005| , Komis 20 - 7 50 keV (|Frederiks fc| 



et al. 20041, and SPI-ACS > 80 keV (Minaev et al. 2009 



20101. Indeed, the T90 parameter of each short burst to be 
summed is less than ~ 2 s. The intensity of the EE seems 
to depend on the energy band being smaller with increasing 



2 M.V.Barkov & A.S. Pozanenko 



100 

Seconds after peak 



Figure 1. Light curve in 50 - 300 keV energy band for 100 short 
(< 1 s), summed, background subtracted, BATSE bursts after 
peak alignment. The time interval corresponding to a peak is not 
shown, adapted from l|Connaughton[2002 ) . 
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Figure 3. A scheme of a two component jet model. Observer will 
register different fluence of short episode of initial emission and 
extended emission depending on the angle of view against axis of 
the jets. 




Figure 2. Light curve of GRB060614 in 15 - 
as recorded with BAT/Swift, time bin is 64 



150 keV energy band 
ms. 



the energy band, and fluence ratio of the EE episode and 
the IPC may be less than 1/100. 

The EE is softer than the IPC, however it is not clear 
yet if this EE is a rising X-ray afterglow or manifestation of 
prolonged activity of the central engine of the burst source. 
Also, based only on the statistical investigation of a large 
amount of short GRBs, one cannot say for sure that each 
particular burst has extended emission. 

Extended emission has actually been observed in indi- 
vidual light curves of some SGRB, confirmed with KONUS 



( Frederiks fc et al.|2004[), HETE -2 ( |Villasenor fc et al.|2005[), 
Swift ( jBarthelmy fc et al.|2005|>, BATSE (jNorris fe Bonnell| 



2006] ), SPI-ACS/INTEGRAL (Minaev et al.||2010[ ). Despite 
the Tgo of some bursts being less than 2 seconds, those the 
EE are significantly detected after re-binning original light 
curve onto larger time scales. In all cases above, the EE 
has a peak flux and fluence much smaller, than analogues 
parameters of the IPC. 

Among the long duration BATSE bursts, several GRBs 



with phenomenological features resembling short bursts 
were found, where IPC does not possess any spectral lag and 
there is a tail of softer emission with durations up to 100 s 
(Norris & Bonnell 2006). The visually estimated amount of 
this burst type is about 2% of the total amount of BATSE 
bursts. While the peak flux of the EE is still 10-30 times 
less than the peak flux of the IPC, the fluence of the EE is 
compatible with the fluence of the IPC. 

Finally, the ultimate example of an apparently long 
burst (Tgo = 102 s, see Figure [2| with all signa tures of short 
bursts is Swift GRB 060614 ( |Gehrels fc etH|2006| |. In ad- 
dition to the absence of spectral lag in the whole light curve 
(both in the IPC and the EE), the host galaxy was identified 
at z=0.125 and deep and temporally dense follow-up optical 
observation did not reveal signature of a supernova, plac- 
ing the brightness of any possible supernova associated with 
GRB 060614 unreasonably faint M R > -13.6 dFynbo & et 



al.|2006|[Gal-Yam fc et al.|2006| [Delia Valle et al.|2006 

It is evident that the intense EE of SGRBs was observed 
in a few cases, whilst the less intense EE can be observed 
more frequently. Using BAT/Swift, which is more sensitive 
to the EE of short bursts, it was shown that among Swift 
SGRB s that ~ 25% of short bursts have EE ( |Norris et aT] 
20111. It is also evident that weak EE are observed in the 



majority of short bursts because sum of the light curves 
possess EE, whilst in individual light curves, the EE is below 
the detection level. 

Thus one could suggest that for short duration bursts 
(or the phenomenon responsible for what we mean as short 
duration bursts) have a distinctive feature, such as extended 
emission with very broad dynamic range of flux and fluence 
of the EE. We suggest that EE is an inherent property of 
short bursts, and an observable property of the EE can be 
explained by the different angular position of the observer 
with respect to the axis of the coaxial jets (see Figure |3|. 
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2 THE MODEL 

In our model, we assume the merging of two NSs or a 
NS+BH system which can lead to formation of a fast- 
rotating black-hole (a — J/M BH , where a is dimension- 
less BH spin parameter) . Simulations of the merging of the 
BH+BH or BH+NS give values of final a as 0.3 < a < 0.65 
( Baker et al.|2008||Ruffert fc Janka|2010| ), and in our calcu- 
lations, we will assume a = 0.5, the formation of relatively 
massive Md ~ 0.001 — 0.2Mq, and compact Rd ~ 10 7 cm 
accretion disk (Lee & Ramirez- Ruiz 20071. Such a disk is 



optically thick and photon cooling is not efficient and in the 
case of ineffective neutrino cooling, such a disk also becomes 
geometrically thick. 

To explain the extended emission of SGRBs, we sug- 
gest a two-component model with neutrino heating ( jWoosley] 
1993 1 and an electromagnetic Blandford-Znajek (BZ) mech- 
anism (iBlandford & Znajek||1977| |Lee et al.||2000| |Mizuno| 



|et al.|2d04||Barkov fc Komissarov|2 008 Rcz zollaet al.|2011| ) 
A short main episode (IPC) of SGRB is the result of a fast 
accretion period (Popham et al. 19991, which launches a 



neutrino powered jet. After a few seconds, neutrino heating 
becomes ineffective, however, the lower accretion rate can 
keep the central engine active for a longer time due to the 
BZ mechanism ( |Barkov fc Ko missarov 2010). The duration 
of the BZ powered jet will depend on the mass in accretion 
disk. 

The accretion time of main mass of the compact disk 
with radius Rd ~ 10 7 cm is short 0.1 — 1 s (van Putten & Os- 
triker||2001[ ). Following [Metzger etH] ( |2008[ ), the accretion 
rate can be estimated as: 



M d « fM d /t m . 



(1) 



where t V i BC = R\jv is viscosity time scale and v is the vis- 
cosity, and the factor / « 1.6. For the viscosity we have 



used an a-prescription (Sha kura||1972| [Shakura fc Sunyaev 
19731, v — ac s H, where c. 



is the isothermal 
sound speed and H is half-thickness of the disc. The initial 
viscosity time scale is: 



0.02aliM- 5 1/2 ^/ 2 



\R d 



(2) 



where M . 5 = M/10°' 5 Af© is mass of the BH, a_i = a/0.1 
is the standard dimensionless viscosity, Rd,7 = Rd/W 7 cm 
is the radius of the accretion disk. 

The accretion disk will be thick and radiatively inef- 
ficient for neutrino cooling if it is larger Rd > 300R g or 
t > Uhick ~ 2q_i s (Chen & Beloborodov 20071, where 
R g = GMbh/c 2 w 5 x lCMo.s cm. For the calculation of 
the extended emission phase, we can use such an approxima- 
tion, whilst t > t t hick, the radioactively inefficient advective 
flux approximation is applicable. We follow the work |Bland-| 
ford fc Begelman] ( |1999[ ), and assume that only a fraction 
~ (R/Rd) p of available material is accreted onto the BH. 
The rest of the mass will be lost to an outflow, and the pa- 
rameter, 'p', can vary from (no wind) up to 1 (powerful 
wind). The self-similar solution for accretion rate was ob- 



tained by Metzger et al. (20081. This solution gives us an 



accretion rate: 



M in « 1.6 



M d .o I R 



tvisc,0 \ Rd,0 



1 + 4.8(1 - C) 



and disk radius 



tvisc,0 



-1 -[l+3(l+2p/3)(l-C)]/[3(l-C)] 



rd « Rd,o 



1+4.8(1 - C) 



t 



tvisc.O 



1 2/3 



(3) 



(4) 



where Md,o, Rd.o and i„i SC ,o are initial mass, radius and 
viscous time of the disk, and C = 2p/2p + 1. The R ma is 



radius of the marginally stable orbit (MSO) (Bardeen et al 



19721 



For the accretion of a thick disk, the maximal luminos- 
ity due to BZ mechanism can be estimated following the 



papers (Komissarov & Barkov 2010 Barkov 20101 where 



used effect of magnetic field grows and a strong large 
scale magnetic field formation in the vicinity of the 



black hole horizon dBisnovatyi-Kogan &j Ruzmaikin 
[19741 [19761 ) ■ 

The pressure of the magnetic field can be a frac- 
tion, 1//3 (where /3 = &itP a /B 2 ), of the gas pressure 
in the disk at MSO, from MSO magnetic field ac- 
cretes to the BH horizon. In such a way, the luminosity 
of the BZ mechanism becomes a weak function of BH spin 
parameter, a, (if 0.5 < a < 1) and can be estimated as: 

0-05 • 2^ 1n 48 -lo-ljV -1 /r\ 

i ^-M in c w 10 a_i/?i M in - 5 erg s , (5) 

a-ipi 



L 



where ft = /3/10. 

The main source of neutrino heating is the neutrino 
annihilation reaction vv — > e + e~. The heating rate can be 
described by ( jZalamea fc Bcloboro dov|2010[ ): 



^« 3 X 10 50 X- 4 - 7 M-^ 2 5 M]f ergs s" 1 , 



(6) 



where X = R ms /4R g . This formula is valid when the ac- 
cretion rate is higher than ~ O.O^a 5 /^ Mq s _1 (Chen & Be 
|loborodov|2007| ) , this critical value of the accretion rate is a 
function of a. As the accretion rate becomes lower, the effi- 
ciency of neutrino heating drops dramatically and becomes 
negligible. 

In our model, the distribution of SGRBs in the inten- 
sity of the extended emission is a selection effect of different 
angular position of the observer (see Figure [3]) with respect 
to the axis of coaxial jets and a dispersion of opening angles 
of the neutrino and BZ jets. The neutrino powered jet has 
the opening angle 6 V „ ~ 0.1 ( |Aloy et"aT . 2005, Harikae et al. 



20101 and can be significantly wider than opening angle of 



1/r (Komis- 



the jet powered by the BZ mechanism 8bz 
|sarov et aL]|2009[ ). 

One can estimate the ratio of jets opening angles us- 
ing GRB 060614. As discussed above, this unique burst has 
brightest EE, i.e. peak fluxes nearly the same in the IPC 
and the EE (Figure^. One can suggest that GRB 060614 
was observed close to the axis of the coaxial jets. Then we 
can write following relation: 



Lb 

L, 



BZ ( 9 v p \ ' 
vv \9bz J 



(7) 



where F\, max and F2, m ax are peak fluxes of the IPC and EE 
episodes. Using equations 1 5|6 1, we obtain 9bz/6vv ~ 0.1. 
Otherwise, one can estimate the same parameter of 9b z /Ovv 
using fluence ratio of the IPC with duration, t vv , of about 
5 s and fluence Si = (16.9 ± 0.2) x 10" 6 erg cm" 2 , and the 
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Figure 4. First episode of short emission and the extended emis- 
sion is observable if initial parameters 'p' and q is in between 
the lines. We take the initial mass of the accretion disk is equal 
to M(j o = O.IMq, a = 0.5 and duration of the IPC and the ex- 
tended emission equals t v p = 1 s, and tsz = 100 s. The area to 
the right of thick line corresponds to Lbz > 3 x 10~ 4 L v p. The 
area to the left of thin line represents the initial accretion rate 
M ln > O.O5M s" 1 



EE episode with duration tsz of about 100 s and fluence 
£h = (3.3 ± 0.1) x 10" 6 erg cm" 2 (see 
( 20061 )): 



Gehrels & et al. 



Lbz 



( Qua y 

\0bz) 



tBZ 



S2 

Si 



(8) 



which leads to the comparable value of 9bz/0 vv ~ 0.1. 

Frequency detection of the EE in the individual light 
curves of short GRB is roughly proportional to the ratio 
of solid angles of the two jets. The frequency ratio can be 
estim ated from 2% of all burst of BATSE ( |Norris fc Bonnell| 
2006[), and up to 25% of short bursts registered by Swift 



iNorris et al 

^Bzj^vv = I 



20111. Hence we have 



t>BZ 



e ui } — h ~=~ 3 and 
15 4- 0.5. One can assume that the opening 
angle ratio of Blandford-Znajek and neutrino powered jets 
is between 0.5 and 0.1. (In above estimations, we ignore any 
possible angular distribution of energy releases and F-factors 
of emitting particles in the jets. We also implicitly assume in 
eqs.(7[8l the conversion factors to the gamma-ray emission 
are equal in both jets.) 



3 DISCUSSION 

The gamma-factors of the two jets can be also different, with 
Y vv ~ 500 predicted from modelling with relativistic hydro- 

12005b. The inter- 



dynamic codes for outflows (Aloy et al 



pretation of observations with a simplified one-zone model 
suggests extremely high Lorentz factor of V > 1200 for the 
IPC of the short GRB 090 510 ( |Ackermann fc et al.||2010| |; 



see also Corsi et al. (20101. Nevertheless, two-zone models 



of GeV photon production gives a significantly lower limita- 
tion on the Lorentz factor of the jets of T > 200 — 400 ( Zou 
et al.|2011 1. 

Ybz can be estimated from the value of open- 
ing angle 9bz ~ 1/Ybz which is specific for BZ-jet 
( |Komissarov et al.||2009[ ), and ratio B z/0 v p < 0.5. If 
Qui, ~ 0.1 we can obtain the lower limit of Ybz > 20, 



which is compatible with numerical calculations of 
BZ-jet parameters at the time the jet is started. In 
the course of BZ-jet development, the Fbz increases 
with time ( jKomissarov et al.|2 009). Then the BZ jet 
leaves the zone of collimating wind, and the outflow 
can get additional acceleration due to the propa- 
gation of a rarefaction wave across the entire jet 



(Tchekhovskoy et al. 2010; Komissarov et al. 2010). 



This additional boosting will not change significantly 
the BZ-jet opening angle while the Tsz could reach 

100 - 300. 

The spectra of the IPC of SGRBs are always harder 
than the spectra of the extended emission and it could be 



explained by Y vv > Ybz (see however |Barraud et al 

(2005)). Another issue is an absence of the spectral lag 
The spectral lag is unavoidable if the emitting particles are 
moving radially (Fenimore et al. 1996 Sari & Piranl 1997; 



Norris 20021. The absence of the spectral lag in the IPC 



may be explained by large value of Y (Norris & Bonnell 



|2006[ ), although, it is not clear whether the lag is present in 
the EE. In a few cases, where the lag could be accurately 



measured, it is negligible within statistical errors (Norris 



& Bonnell 2006 Gehrels & et al. 20061. In general, more 



investigations of phenomenology and nature of the spectral 
lag are necessary. 

The IPC of GRB 060614 falls far off (IGehrels & et al. 



2006 1 the Amati relationship for long GRBs ( Amati & et 
al.|2002 Amati|2010 1. However, the E p , E iso parameters of 



the whole burst are consistent, within 2 sigma level, with 



the Amati relationship for long GRBs (Mangano & et al. 



20071. One would observe the GRB 060614 off-axis of the 



BZ-jet, such could not detect the EE, but only IPC. In this 
case, eventually GRB 060614 would be an outlier of the Am- 
ati relationship as we see it for other short duration bursts 
( Amati|2010| . 

Finally let us estimate physical parameters when the EE 
is observable. Based on the equation (JsJ) and typical values 
obtained for extended emission registered for sum of the 
short duration bursts, i.e. fluence ratio S2/S1 ~ 1, t vu ~ 1 s, 
tBZ ~ 100 s and mean value of 9bz/0 vv ~ 0.2, we can infer 
that the extended emission can be detected if Lbz/L^h > 
3 x 10" 4 . 

Using eq.|3| at time t = 0, we can draw in the Figure[4|a 
(thin) line which shows the limit max(Mi n ) = 0.05M Q s~ . 
Thick line represents the limit Lbz/L^d = 3 x 10 -4 , and 
the area to the right of this line implies that there is enough 
energy in the BZ jet to be detected. The region between 
the lines provides a region in which both IPC and extended 
emission could be detected. Indeed in the region on the left 
of the thick line (Lbz/L vu < 3 x 10 -4 ), the accretion rate 
is ether too slow to start BZ jet or to provide enough energy 
in the EE. 

On the other hand, in the region to the right of thin 
line, the initial short period of a fast accretion rate provid- 
ing a neutrino powered jet is absent. However, the accretion 
rate can be sufficient to start the BZ jet. Therefore one can 
expect the existence of rare type of long duration bursts 
which look like only the extended emission of such bursts as 
GRB 060614. The observable features of the burst are the 
absence of a supernova, negligible spectral lag, and a pos- 
sible precursor corresponding to a failed neutrino jet. It is 
difficult to estimate a frequency of the occurrence of such 
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bursts as it will depend on the initial accretion rate dis- 
tribution. However, one can remember that ~ 10% of long 
BATSE bursts (Hak kila fe et al.|2007| > actually possess zero 
spectral lag (see also (Norris 2002)). Indeed, some of them 
are short bursts with observable extended emission, as found 
by jNorris fe Bonnell||2006"l ), whilst some of them could be 
only BZ powered bursts. Recent spectral lag investigation 
of naturally long burst of the BAT/Swift experiment gives 
the same ~ 10% of burst which have negligible lag within 
statistical error ( |Ukwatta et~aI1[20To| . 

The main source of the magnetic field which is neces- 
sary to start the BZ jet can be a seed magnetic field of the 
merging companions. The variability of magnetic flux can be 
a manifestation of a magnetic dynamo in the accretion disk 



and the accretion of magnetic flux onto the BH (Barkov & 
Baushev]|2011| [Rezzolla et al.||2011| ). The accretion of mag- 
netic flux with different polarity leads to the variability of 
magnetic flux on the BH horizon and thus a variability of 
BZ luminosity as well. 

The magnetic field generated due to a dynamo can have 
a pressure, which can be comparable with the pressure of 
matter ( Bisnovatyi-Kogan & Lovelace 200 7] |Rothstein fc| 



Lovelace 2008 Igumenshchev 2008 Lovelace et al. 2009) 



The pressure of the gas in the disk is a function of its thick- 
ness, P oc (H/R) 2 , and hence B oc (H/R), whilst the energy 
release will be Lbz oc B 2 oc (H/R) 2 . On the time scale of 
t < tthick, the magnetic dynamo can generate significantly 
weaker fields than in the thick disk and luminosity of BZ 
is suppressed. In the time scale of t > tthick which may be 
of few seconds the disk becomes thicker, the magnetic dy- 
namo becomes more efficient and provides an increased BZ- 
jet luminosity. The increase of the luminosity may explain 
the hiatus between the IPC and the extended emission ob- 
served in some short burst with EE ( jNorris fc Bonnell|2006| 
Gehrels fc et al.|2006 l. 

Finally, the same mechanism involving two jets can be 
responsible not only for GRB based on NS-NS or NS-BH 
merging. Indeed the merging of BH or NS and white dwarf 



( Fryer et al. 1999 1 or BH and a core of Wolf-Rayet star ( Fryer 
fc Woosley|1998j > can rise a long duration burst and exhibit 
a long X-ray afterglow (Barkov & Komissarov 20101 with 



appropriate time scaling of the main episode powered by a 
neutrino jet and extended emission, powered by BZ jet. The 
time scaling parameter will be dependent on the initial linear 
size (specific angular momentum) of the accretion disk, and 
in the case of a BH and white dwarf merging, may be ~ 10. 
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